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p-HARMONIC COORDINATES FOR HOLDER METRICS 
AND APPLICATIONS 


VESA JULIN, TONY LIIMATAINEN, AND MIKKO SALO 

Abstract. We show that on any Riemannian manifold with Holder 
continuous metric tensor, there exists a p-harmonic coordinate system 
near any point. When p = n this leads to a useful gauge condition for 
regularity results in conformal geometry. As applications, we show that 
any conformal mapping between manifolds having C a metric tensors is 
C 1+ “ regular, and that a manifold with W 1,n fl C a metric tensor and 
with vanishing Weyl tensor is locally conformally flat if n > 4. The 
results extend the works ILS141ILS15I from the case of C 1+a metrics 
to the Holder continuous case. In an appendix, we also develop some 
regularity results for overdetermined elliptic systems in divergence form. 


1. Introduction 

Various regularity results in Riemannian geometry can be achieved via 
the use of harmonic coordinates. If ( M,g ) is a Riemannian manifold, a local 
coordinate system x = (x 1 ,..., x n ) is called harmonic if 

A g xi = 0, 1 < j < n, 

where A g is the Laplace-Beltrami operator. The notion of harmonic coordi¬ 
nates goes back to Einstein |Eil6j . and these coordinates have the property 
that the Ricci and Riemann curvature tensors become elliptic operators 
when written in harmonic coordinates DK81 . The Ricci and Riemann ten¬ 
sors are invariant under diffeomorphisms, and fixing a harmonic coordinate 
system may be thought of as a gauge condition which makes these opera¬ 
tors elliptic. Harmonic coordinates are useful when studying regularity of 
isometries or flatness of low regularity metrics | Ta06] . Also, the DeTurck 
trick for local wellposedness of the Ricci flow ('KOI is related to harmonic 
coordinates, see (rmrrj . 

Harmonic coordinates are well suited to studying regularity results in Rie¬ 
mannian geometry, due to the fact that harmonic functions are preserved by 
isometries. In this paper we are interested in regularity results in conformal 
geometry. For this purpose it is useful to have a gauge condition which is 
preserved by conformal transformations. Such a condition was introduced 
in the recent work }LS14| . which established the existence of p-harmonic 
coordinates on any smooth manifold with C s , s > 1, Riemannian metric. 
If p = n, the corresponding n-harmonic gauge condition is conformally in¬ 
variant, and II,S1 11 applied n-harmonic coordinates to give a new proof of 
certain regularity results for conformal mappings. The subsequent article 
[LS15] showed that n-harmonic coordinates provide a gauge condition for 
ellipticity of conformal curvature tensors, and characterized local conformal 
flatness of manifolds of dimension > 4 having C s , s > 1, metric tensor. 
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The aim of the present article is to extend the results of [LS14 ] to C s , 
s > 0, metric tensors and the local conformal flatness characterization of 
[LS15] to W l ' n D C s , s > 0, metric tensors. Let us state the desired results 
in more detail. 

In this paper M will be a smooth (= C°°) manifold of dimension n > 2, 
and g will be a Riemannian metric on M with finite regularity. If 1 < p < oo, 
the p-Laplace operator on ( M,g ) is the operator u i-A — 5{\du\ p ~ 2 du) where 
d is the exterior derivative, S is the codifferential and \du\ is the g -norm of 
the cotangent vector du. This operator is given in local coordinates by 

- 5(| dur 2 du) = Igr^djQg^^ig^daudbU^dku) (1.1) 

where (gjk) is the metric in local coordinates, (g jk ) is the inverse matrix 
of ( gjk ), and \g\ = det (gjk)- Throughout the paper we use the Einstein 
summation convention, where a repeated index in upper and lower position 
is summed from 1 to n. If p = 2,0 is the usual Laplace-Beltrami operator, 
whereas for p / 2 the p-Laplace operator is a quasilinear degenerate elliptic 
operator. We refer to (HKM93I iLinOG) for more details on p-Laplace and 
more general „4-harmonic operators. 

We will state our results in terms of C s spaces for s > 0. If s is a 
nonnegative integer, we understand that C s (M n ) is the space of functions in 
M n whose derivatives up to order s are bounded and continuous. If s = k + a 
where k > 0 is an integer and 0 < a < 1, then C s = C k+a is the usual Hblder 
space with norm 


ll/llc fe +' 


V ll^/llioopH-n) + SU P 

l7l<fc h\=k 


\d^f(x) -d^f{y)\ 
\x — y\ a 


We will also consider the Zygmund spaces C*(R n ) for s G R, which are 
defined via the norm 


II/IIci = sup 2 3S \\lpj(D)f || L oo (K n) 
j>0 

where (V’j(O)^o is a Littlewood-Paley partition of unity. One has 
C s = C®, s > 0 not an integer, 

and C k C C k if k is a nonnegative integer. Most often we will use local 
versions of the above spaces, and in local coordinate statements will write 
C s instead of Cj^ c etc. See |Ta96] for more details on the above spaces. 

The first main theorem states that whenever (M, g) is a Riemannian man¬ 
ifold with C s metric tensor where s > 0, then near any point of M there 
exist local coordinates all of whose coordinate functions are p-harmonic. 
This extends a result from (LS14) from the case s > 1 to s > 0. 


Theorem 1.1 (p-harmonic coordinates). Let ( M,g ) be a Riemannian man¬ 
ifold whose metric is of class C s , s > 0, in a local coordinate chart about 
some point xq G M. Let also 1 < p < oo. There exists a local coordinate 
chart near xo whose coordinate functions are p-harmonic and have C| +1 
regularity, and given any e > 0, one can arrange so that metric satisfies 
\gjk{x o) — Sjk\ < £ for j, k = 1,..., n. Moreover, all p-harmonic coordinates 
near xq have C* +1 regularity. 
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The next result states that if 0 is a 1-quasiregular map between two 
manifolds with C s metric tensors where s > 0 (in particular i f> could be a C 1 
conformal map), then </> is in fact a C s+l local conformal diffeomorphism. 
This extends the results of |Iw82l l.SlT from the case s > 1 to s > 0. The 
earlier result |Le76| considered C°° Riemannian metrics, and [Ta06l ICL15J 
give analogous results for Riemannian and subRiemannian isometries. 

Theorem 1.2 (Regularity of conformal maps). Let (M,g) and ( N,h ) be 
Riemannian manifolds, n > 3, where g,h E C s for some s > 0, s / 1. Let 
4> : M —> N be a non-constant mapping. Then the following are equivalent: 


0 is a Riemannian 1-quasiregular mapping, (1.2) 

<f> is locally bi-Lipschitz and ip*h = eg a.e., (1.3) 

<j> is a local C 1 diffeomorphism and (f>*h = eg, (1.4) 

(f> is a local C s+1 diffeomorphism and <j)*h = eg. (1.5) 


The third result states that the vanishing of the Weyl tensor on a manifold 
with W 1,n C\C a metric tensor is equivalent with local conformal flatness. The 
Weyl tensor W ( g) of (M, g) is the 4-tensor 

W^abcd Robed T Pac9bd Pbc9ad T Pbd9ac Pad9bc 

where R a bcd = g(R(d a ,db)d c ,d c i ) is the Riemann curvature tensor, and P a b 
is the Schouten tensor given by 

Pab 7T I Pab 7TT T\9ab ) 

n — 2 \ 2 (re — 1) ) 

where Rb c = g ad R a bcd is the Ricci tensor and R = g rs R rs is the scalar 
curvature. If n = 3 the Weyl tensor is always zero, and if n > 4 the Weyl 
tensor is conformally invariant in the sense that 

W(cg) = cW(g) 

if c is a smooth positive function. It is a classical result that for C 3 metrics, 
the vanishing of the Weyl tensor is equivalent with the manifold being locally 
conformally flat in dimensions n > 4 (see for instance [Au88l Chapter 4]). 
The question of having a low regularity counterpart for this result has been 
raised in [IMP 11 Section 2.7] and |Mal4l Section 1], Such a result was proved 
for C s metrics with s > 1 in [LS15 ]. Here we improve it to W 1 ' n C\C s metrics 
with s > 0, where W k,p is the usual (local) Sobolev space of functions whose 
derivatives up to order k are (locally) in L p . 

Theorem 1.3 (Local conformal flatness). Let M be a smooth n-dimensional 
manifold with n > 4, and let g E W 1,n n C a in some local coordinates near 
xq E M where 0 < a < 1. If W (g) = 0 near xo, then g a b = c5 a b for some 
positive W 1,n (~l C a function c in some ra-harmonic coordinates near xq. 

In three dimensions, local conformal flatness of sufficiently regular metrics 
is characterized by the vanishing of the Cotton tensor. An analogue of 
Theorem P for n = 3 was proved in [LS15| for C s metrics with s > 2. It 
is likely that the regularity assumption could be improved slightly, but we 
will not consider this here. 
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We also remark that there are alternative gauge conditions for ellipticity 
of conformal curvature tensors, based on solving the Yamabe problem lo¬ 
cally to achieve constant scalar curvature and using harmonic coordinates. 
This procedure may require more regularity than the method of n-harmonic 
coordinates presented in this paper, see [LS15j for further discussion. 

This paper is organized as follows. Section |T] is the introduction. In 
Section [2] we prove Theorem ll.il and more generally we show the existence of 
^4-harmonic coordinates for any A satisfying certain structural assumptions 
including Holder continuity in x. The main point in extending the result of 
}LS14] from s > 1 to s > 0 is that instead of working with strong solutions 
and nondivergence form equations, we keep the equations in divergence form 
and apply Campanato and Schauder estimates in order to show regularity 
of weak solutions. In Section [2] we also give an additional W 2,q fl C l+a 
regularity result for the p-harmonic coordinates if g is W 1 ’ 9 fl C a where 
0 < a < 1 and q> 2. 

Section [3] considers the regularity of conformal mappings and establishes 
Theorem ll.2l Given Theorem ll.il the argument is analogous to [LS14] and is 
based on representing the conformal map locally in n-harmonic coordinates 
in the target manifold. We also supply some details related to integer values 
of s that were missing in [LS14j . 

In Section [4] we discuss the Weyl tensor for low regularity metrics and 
prove Theorem ll.3l This uses the fact proved in [LS15j that the Weyl tensor 
becomes elliptic in n-harmonic coordinates after a conformal normalization. 
The vanishing of the Weyl tensor can be interpreted as a quasilinear overde- 
termined elliptic system involving a nonlinearity with quadratic growth, and 
to deal with this we will require regularity results for overdetermined elliptic 
systems. 

There are many regularity results for square elliptic systems available in 
the literature, see for instance |Gi83l CM 12 . In the case of sufficiently 
smooth coefficients, a linear overdetermined elliptic system Pu = / can 
be reduced to the square system P*Pu = P*f, and regularity results will 
follow from the corresponding results for square systems. Some results for 
overdetermined systems with nonsmooth coefficients (mostly in nondiver¬ 
gence form) are given in [Ta961 ITaflOl ILS15] . However, we need to deal with 
certain systems in divergence form which do not seem to be covered by these 
results. For this reason, in Appendix [A] we establish the basic Calderon- 
Zygmund and Schauder regularity results for linear overdetermined elliptic 
systems in divergence form. Complete proofs for these results are included, 
based on freezing coefficients and Fourier arguments. 
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2. p-HARMONIC COORDINATES 

In this section we will prove Theorem 11.11 Since this theorem is a local 
statement, it is sufficient to consider a corresponding result in M n . Locally, 
p-harmonic functions in (M, g ) are solutions of the .4-harmonic equation 


div A(x, Vu) = 0 (2-1) 

in some open subset 12 C M n where 

A j (x,0 = \g(x)\^g^(x)(g a \x)U b ) P -^^ ( 2 . 2 ) 

The next result shows the existence of .4-harmonic coordinates for any A 
which satisfies the following structural assumptions for some M, 6 > 0: 

A(x,tQ = 2 p_1 .4(x,£) for t> 0, (2.3) 

\d^A k (x,0\<M\cr 2 , (2.4) 

\A{x,£) - A(y,0\ < M\x - j/|“|C| p_1 J (2.5) 

{A{x, 0 - A(x, 0) • (£ - C) > <KI£I + ICir 2 |^ - Cl 2 - (2-6) 


Theorem 2.1 (.4-harmonic coordinates). Let 12 be an open set in M n and 
let 1 < p < oo and 0 < a < 1. Suppose A : 11 x M 71 —>• is a vector- 

held which satisfies (12.3D - (12.61) . Given any point xq £ 12, there is an C 1+a 
diffeomorphism U from some neighborhood of xq onto an open set in M n 
such that all coordinate functions of U are .4-harmonic. Moreover, any C l 
diffeomorphism whose coordinate functions are ^4-harmonic has this regu¬ 
larity. Given any invertible matrix S and any e > 0, there is U dehned near 
xq such that 

\\DU(x 0 ) - S\\ < e. 

We remark that (12.3D (12.6D imply the conditions (for x £ 12 and £, h £ M n ) 

•A(®,0-£> W, 

d^{xA)h j h k >8\ir 2 \h\\ 

which appear in the C 1+/3 regularity results for the .4-harmonic equation 
(ED |Ma86j . 

We follow the argument of [LS14j except for one part of the proof. In 
mM, we used the fact if A is C s in the x variable where s > 1, then an A- 
harmonic function whose gradient is nonzero must be C' s+1 . The assumption 
s > 1 was used to first show that the .4-harmonic function is IT 2 ’ 2 by a dif¬ 
ference quotient argument, and then the .4-harmonic equation was brought 
into nondivergence form where Schauder estimates could be applied. In this 
paper, we will work under the assumption s > 0 and apply Campanato type 
estimates and a frozen coefficient argument to the divergence form equation 
instead. 

We recall the following classical result by Campanato. To that aim we 
denote the mean value of it in a ball B(xo,r) by u xo , r . The proof can be 
found in [Git)3] . 
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Proposition 2.2. Suppose that u £ L 2 (Q ) satisfies 

I \u — u X:r | 2 < Cr n+ 27 for every B(x,r) C D, r < ro, 

J B(x,r) 

for some 7 € (0,1) and ro > 0. Then u £ C* 7 (D). 

We will also need the following standard algebraic lemma (see [CM 
Lemma 7.3]). 

Lemma 2.3. Let 0 be a nonnegative and nondecreasing function on [0, ro] 
such that 

4>{p) < A (-'j <j)(r) + Br 13 for every 0 < p < r < ro, 

where A, B, a, (3 are nonnegative constants such that (3 < a. Then for every 
7 £ (/?, a) there exists a constant C such that it holds 


<£(p) < c ((~) 7 < M r ) + b p p 


for every 0 < p < r < ro- 


The next result concerns C 1+/3 regularity of .4-harmonic functions, and 
requires the assumption that A is C a with respect to x for some a > 0 . 
We recall that even in Euclidean space, solutions of the p-Laplace equation 
always have C 1+/3 regularity for some (3 > 0, but f3 could be very small and, 
in particular, solutions are in general not C 2 regular. 

We write B r = B(0,r) for the open ball of radius r > 0 centered at the 
origin in M n , and [f] cP = sup x ^ y ■ 

Proposition 2.4. Suppose u £ W l,p (B(xo,r)) is a weak solution of 

div.A(z, Vu) = 0 

in B(xo,r) and A satisfies (12.3D112.6p . Then there are constants C > 0 and 
f3 £ (0,1) only depending on n, p, a, 5 and M such that u is C 1+/3 regular 
in B(x 0 , r) and 

l|VM||L°°(B(x 0 ,r/ 2 )) +^N C ^ B ( 10ir / 2 )) < Cr n,l ’\\Vu\\ L p( B ( X0>r )). 

Proof. We may assume that xq = 0 and u solves div.4(x,Vr() = 0 in B r . 
Define scaled coordinates x = x/r and a scaled function u(x) = r^ 1 u(rx). 
Then u £ W 1,P (B\) is a solution of 

div A r {x. Vu(x)) = 0 in B\ 

where A r (x,f) = A(rx,f). Now A r satisfies ()2.3D (12.611 with the same con¬ 
stants 5 and M (this uses that r < 1). 

We define 

u - n 0 1 

v = —---- 2 -. 

\\u ~ 7to,i||w 1 .P(B 1 ) 

Then v £ W l,p {B\) is a solution of div.4 r (x, Vv(x)) = 0 in B\ with 
ll^llw^plSi) = 1- If follows from |Ma 86 j (see also [ Se(i 1 . [ DiB83] ) that 
there are C > 0 and (3 > 0, only depending on n,p, a , 5 and M, such that 

IMIc^l b 1/2 ) < c. 
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In particular, this implies that 

l|V«llL~ ( B 1/a ) + \^ V \cP{B 1/2 ) < C- 

The Poincare inequality yields ||{t — uo,i \\w 1 ^(B 1 ) — C\\V , u\\ L p(b 1 )- Conse¬ 
quently 

l|Vu|| L< x»( Bl/2 ) + [Vu] cHBl/2) < C\\Vu\\ LP ( Bl y 

Undoing the scaling proves the result. □ 


The following lemma will be useful for the Schauder estimate that we 
use later. We will apply it in the case where either the solution of the A- 
harmonic equation or the solution of the equation with frozen coefficient has 
nonvanishing gradient, and L 2 estimates will be natural in this context. 

Lemma 2.5. Suppose A satisfies (J2.3IMJ2.6I) and u\,uo € W 1,p (B(xq, r)), 
r < 1 , are such that u\ — uq £ Wq' p (B(xq, r)) and u\ is a weak solution of 

div*4,(x, Vrti) = 0 

and uo is a weak solution of the equation with coefficients frozen at xq, 

div A(xo, V'Uq) = 0 . 

If either u\ or uq satisfies the condition 


c 0 < |V«i| < 1/co in B(x 0 ,r) (2.7) 

for some Co > 0 , then we have the two estimates 

||Vui - Vu 0 \\l 2 ( B ( Xo ,t/ 2 )) < Cri +a , ( 2 . 8 ) 

2ol/3 

llViii - Vu 0 ||ioo( B(a . 0 ,r/4)) < (2.9) 

The constants C and /3 depend on co, n, p, a, 5 and M. 

Proof. Without loss of generality we may assume that xo = 0. We prove 
the result by assuming that uo satisfies m since in the other case the 
argument is similar. Write A\ = A and Ao {x,0 = -4(^0, 0 - 

Let us first prove (12.81) . We use the assumptions on A and the fact that 
we may use u\ — uq as a test function in the equations for u\ and ilq to 
conclude that 


/:= / (|Vui| + |Vti 0 |) p_2 |Vui — Vu 0 | 2 dx 

{A\(x, Vui) — Ai(x, Vuo)) • (Vui — Vuo) dx 
-iW 1 / (Ai(x, Vuo) — Ao(x, Vuo)) ■ (Vui — Vuo) dx 

J B r 


' B r 

< 5' 1 


' B r 


< 8 


-l 


' B r 


(|V«i| + |Vu 0 |) 2 p \Ai(x, Vu 0 ) - A 0 (x,Vu 0 )f 2 dx 


1/2 


I 1 /2 





p-HARMONIC COORDINATES FOR HOLDER METRICS 

Thus by (12.51) 

I < 5~ 2 
M 2 


< - r 

~5 2 


B r 

2a 


(|V«i| + |Vn 0 |) 2 - p |^i(x, Vn 0 ) - A 0 (x, Vn 0 )| 2 

/ (|Viti| + |Vno|) 2_p |V-u 0 | 2p_2 . (2.10) 

J B r 

Assume first that p > 2. Then (12.71) and (12.101) imply 

[ |Vui - Vu 0 | 2 < Cr 2a [ |Vn 0 r 
J B r /2 ** 

which is (12.8|) . 

Let us assume 1 < p < 2. We estimate the right hand side of (12.101) by 
Holder inequality with exponents q = and q' = 2^2 aR d obtain 


' B r 


< 


(|Vui| + |Vn 0 |) 2 " p |Vu 0 | 2p - 2 


2-P 

V 


B r 

2 —pi 


|VuoP 


2p—2 

P 


( 2 . 11 ) 


(|Vm| + |Vu 0 |) p 

3r 

< (l|V«i|| LP(Sr ) + \\Vu 0 \\ LP{Br )) 2 -^u 0 f L p - 2 Br) 

Moreover, since f B A\(x, Vu 1 ) • V(iti — uq) dx = 0, we have 

5 / \Vui\ p dx < / A\(x, Vrti) • Vui dx = / A\{x,\7u\) ■ \7 uq dx 
J B r J B r J B r 

<M j iVtiiMVtiol dx < MllVuif-^jVuohptBr) 


Hence 


M, 


l|Vni|| LP( Br ) < -rl|V«o|| LP(B r )- 


( 2 . 12 ) 


Therefore the the right hand side of (12.101) can be estimated by (12.111) and 

(Era 


2 a 


[ (|Vni| + |V« 0 |) 2 “ p |Vnol 

J B r 


2p ~ 2 < Cr 


2a 


' B r 


|Vn 0 | ? 


(2.13) 


To estimate the left hand side of (12. 10[) we use Proposition 12.41 and (12.121) 
show that on B r /2 

|Vuo| + |V«!| < 1/co + Cr- n l p \\y Ul \\ LP{Br) 

<l/c 0 + Cr- n / p \\X7u 0 \\ LP(Br) 

< C. 

Therefore (12.101) and (12.131) yield 

|Vtti — Vn 0 | 2 < Cr 2a f |Vn 0 r 
) J B r 


' B, 


/ 2 


and we again have (J2.8I) . 
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It remains to prove (12.91) assuming (12.71) for uq. We first invoke Proposi¬ 
tion [23] which shows that 

AVuj] cnBr/2) < Cr~ n / p 11 Vuj\\ lp{b t )i 3 = 0,1. 

Using () 2.121) and the fact that | Vito I < 1/co in B r , we obtain 

r P [V Uj \ c P {Br/2) <C, J = 0,1- 

Combining this with (12.81) and increasing the constant slightly gives the two 
estimates 


[^(Vn 1 -Vn 0 )] c / 3(Sr/2) <C4f+^ 

||^(Vm - V« 0 )IU 2 (B r/a ) < 

Interpolating these two estimates by [LS141 Lemma A.l] (with the choices 
II = B r j 2, K = B r / 4 , M = CAz +l3 and / = r^(Vui — Vito)) gives 


r^Vui - V«o)||L»(B r/4 ) 


< C n!Pi/3 (C4t +/3 )^W||r /3 (Vni 


< Cr P+ 


2 cep 
n+2/3 # 


Vn 0 )|| 


2/3 

n+2/3 


L 2 (B r/2 ) 


This proves (12.91) . 


□ 


The next result shows that if A is C a with respect to x, then one can 
upgrade the C 1+ @ regularity of a weak solution u (recall that (3 could be 
very small) to C l+a regularity in any open set where Vit is nonvanishing. 

Proposition 2.6. Suppose it £ W l ' p (B\) is a weak solution of 

div A{x, Vit) = 0 

and A satisfies (|2.3D - (|2.6D . If Vit is nonvanishing in some open set U C B i, 
then u is C l+a in U. 


Proof. Assume that Vit is nonvanishing in U C B\. Our aim is to show 
that u is C 1+a near any xq £ 17. Fix xq £ U and choose R > 0 so that 
B(x o, 8 R) C U. We will show that for some small ro > 0 one has 

[ |V» - (V«) 2 , r | 2 dy < Cr"+ 2 «||V«|| 2 2(B(I0 , 8S)) (2.14) 

J B(z,r) 

whenever B(z,r) C B(xo,R ) and r < ro. This implies it £ C 1+a (B(xo, R)) 
by Proposition 12.21 

By Proposition 12.41 it is C 1+ P regular in B\ and thus there is a constant 
Co > 0 such that 

Co < |Vit| < 1/co in B(xo,8R). (2.15) 

Fix z and r with B(z, r) C B{x o, R), and let v = v r £ W 1)V {B( y z, 8r)) be the 
weak solution of the frozen coefficient equation 

&ivA(z, Vu) = 0 in B(z,8r), v = u on dB(z,8r). 
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By Lemma 12.51 we have (with C independent of z and r) 

l|Vu - Vv|| L 2 {B(2;4r)) < Cri +Q , (2.16) 

l|V« - Vu|| L oo ( B (z, 2r)) < Cr™+2P. (2.17) 

The bounds (12.1511 and ()2.17f) imply that there is ro > 0 depending on co 
such that 

|Vu| > co/2 in B(z,2r) (2.18) 

when r < ro- Now v is the solution of an .4-harmonic equation in B(z,2r) 
with constant coefficients and Vv is nonvanishing by (12.1811 . A standard 
argument (see for instance [LS14L Proposition 2.3]) implies that v is C°° in 
B(z,r). Thus arguing as in |Gi031 Theorem 7.7] we have the estimate 

[ \Vv-(Vv) z J 2 dx<c(^) n+2 f |Vu - (Vu ) 2>r | 2 dx (2.19) 

J B(z,p) Vr/ J B(z,r) 

for every p < r < ro for some small r o- We use (12.161) . (12.191) and the 
estimate f B ^ z ^\f — f z ,p\ 2 < f B (. p ^\f — A] 2 for any A £ M to conclude that 


f | Vu-(Vu) z>p \ 2 dx 
J B(Z'P) 


< 2 


< 


/ |Vu — (Vu ) 2iP | 2 dx + 2 / |V(tt — v)\ 2 dx 

I B(z,p) J B(z,p) 

C(T 2 I |Vv-(Vu) 2 , r | 2 dx + 2 [ \V(u-v)\ 2 dx 

A r' J B(z,r) ’ J B(z,r) 


< Cl 


B(z,r) 

fpy +2 f |Vu — {Wu) Z) r\ 2 dx + C [ |V(u — v)\ 2 dx 

A T 7 J B(z,r) ’ JB(z.r) 


< c ( p -) n+2 


<B(z,r) 


B(z,r) 

\Vu- (Vu) z , r \ 2 dx + Cr n+2a 


for every p < r < ro- The estimate ()2. 141) then follows from Lemma 12.31 □ 


We continue to study the regularity properties of 4-harmonic functions 
and prove a Calderon-Zygmund type result which will be needed later in 
Section [2 We denote the partial derivatives of x 4 4(x,£) by d Xi A(x ,£) 
and of £ i-a A(x,£) by d^A(xA)- Let us assume that in addition to the 
structural assumptions (|2.3I) - (|2.6D the vector field A satisfies 


1 d Xi A(x,£)\ < f{x )|£| p S 

for some / £ L 9 (fl), 

( 2 . 20 ) 

where q > 2 and 



(x,0 ^ %-4(x,£) 

is continuous. 

( 2 . 21 ) 


Proposition 2.7. Suppose u £ W l ' p (B\) is a weak solution of 


div4(x, Vu) = 0 


and A satisfies (12.31) ([ 2 . 6 D and ( 12 . 201 ) . ( 12 . 211 ) . If Vu is nonvanishing in some 
open set U C B±, then u £ W 2 ^(U). 
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Proof. Assume that Vu is nonvanishing in B(xo,2r) C B\. Without loss of 

generality we assume that x$ = 0. By the previous theorem we have that 

u £ C 1,a (B3r). We show first that u £ W 2,2 (B r ). 

2 

Let us fix a direction e* and choose a cutoff function ( £ C^°(B 4 r) such 

that 0 < C, < 1 and £ = 1 in B r . For a small h > 0 denote the difference 
quotient by A hu(x) = \(u(x + hef) — u(x)). We use the test function 
ip = A hu(x)C 2 and integrate against the equations 

—div*4,(x, Vu(x)) = 0 and — div.4(x + hei, Vu(x + kef)) = 0. 


By rearranging terms, this gives that 

(A(x + hei , Vk(x + hei )) — A(x + hei, Vu(x))) ■ (VA hu)( 2 dx 


B 4r 

~T 


< 


/ |*4(x + hei, V'u(x)) — .4(x, Vit(x))| |V Ahu\( 2 dx 

T 

+ 2 / |„4.(x + /re*, Vu(x + /le^)) — A(x + hei, Vu(x)) J | Ahu\(\ V£| dx 

J 5 4t* 

T 

+ 2 / |^4(x + hei, Vu(x)) — „4(x, Vu(x))| |A/,u|C|VC| dx. 


’ B 4r 

~T 


( 2 . 22 ) 


We use the assumption (|2.6I) . the C 1,a regularity of u and the assumption 
|Vn| > c > 0 to estimate the left hand side of ()2.22l) 


' B 4 r 


(A(x + hei, Vu(x + hei)) — A(x + hei, Vrt(x))) • (VA hu)( 2 dx 


> 


ch f |VA hu\ 2 ( 2 dx 

•J B 4v 


(2.23) 


for some constant c > 0. Next we use (12.201) and Young’s inequality to 
estimate the first term on the right hand side of (12.221) 

/ \A{x + hei , Vtt(x)) — A(x, V«(x))| |VA/,n|C 2 dx 

)B 4t 

T 

< — f |VA hii\ 2 () 2 dx + ^- f \A(x + hei, Vu(x)) — A(x, Vit(x))|" dx 
4 JB ir d J B 


5 4 r 

"3" 


> 4r 
"3" 


< 


< 


ch 


[ |VA h u\ 2 C 2 dx + ^ j j — A(x + tei, Vu(x)) dt 

•J B 4 r ** B 4 r J 0 

~T T 

fh 


dx 


C -B f |VA/jti| 2 C 2 dx + C f f \d Xi A(x + tei, Vu(x ))| 2 dxdt 
^ J B 4 r J 0 J B 4 r 

~T T 


< 


ch 


j |VA h u\ 2 ( 2 dx + Ch f f 2 dx. 

B 4 v ** B 

~T 


B 3r 
~T 


(2.24) 
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Similarly we have 

\A(x + hei, Vu(x)) — A(x, Vit(x))| |A/jit|()|V£| dx 


B Ar 

~T 


< 


Ch [ (|Vu| 2 + f 2 ) dx. 

J B 3 r 

~T 


(2.25) 


Moreover by the assumption (12.41) . the C 1,a regularity of u and |Vit| > c > 0 
we have that 

/ | A(x + hei, Vu(x + hei)) ~ A(x + hei, Vu(x))| |Ahu|C|VC| dx 

J B 4r 

cf \Vu(x + hei) - Vu(x)\ \A h u\(\V(\dx ( 2 .26) 

J B 4r 
"3" 

^ [ \VA h u\ 2 ( 2 dx + Ch [ |Vu| 2 dx. 

^ JB A r JB 3 r 


>4r 

"3" 


5 4r 

"T 


>3r 

T 


We divide (12.22(1 by h and deduce from (12.231) . (12.241) . (12.251) and (12.261) 
that 


' B, 


|VA h u| 2 dx < C f (/ 2 + |Vu| 2 ) dx. 

J Bzr 

~T 


By letting h —> 0 we obtain that u G W 2,2 (-B r ). 
We may thus differentiate the equation 


—div„4(x, V«(x)) = 0 

with respect to x* and conclude that v = d Xi u is a solution of 

— div (d^A(x, Vu(x))Vv) = div (d Xi A(x, Vit(x))) . (2.27) 

By the assumption (12.61) the linear operator on the left hand side of (12.271) is 
elliptic. Moreover by ([2.211) x i-A d^A(x, Vtt(x)) is continuous. On the other 
hand the vector field x i-A d Xi A(x,Vu(x)) on the right hand side of (12.271) 
is L q integrable and therefore by the classical Calderon- Zygmund estimate 
(also contained in Proposition IA.H all we have that v G W l,q . This implies 
u G W 2 ' q {B r ). □ 


Proof of Theorem \2.1[ We may assume that xo = 0, B\ C and S is the 
identity matrix (otherwise consider ( Sx) J instead of x 3 below). The A- 
harmonic coordinates will be obtained as the map U = (it 1 ,... ,u n ), where 
each u 3 = ui solves for r small the Dirichlet problem 


divM(x, ViP) = 0 in B r , u j \dR r = x- 7 . 

Let us fix j G {1,..., n}. Then v 3 (x) = x- 7 is a solution of divM(0, Vu- 7 ) = 0 
in B r and 


\Vv J \ = 1 in B r . 
We use Lemma [231 to conclude 


2 ap 


IIVxc 7 - Vu J || L oo (Sr/4) < Cr"+2^ 
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Thus there is r > 0 such that |Vu 3 — Vi4| < ^ in B r / 4: and, in particular, 
ViP is nonvanishing in B r / 4 . Proposition 12.61 shows that u 3 £ C 1+a (B r / 8 ). 

Now each u 3 is C 1+a near 0 and the Jacobian matrix of U = (it 1 ,..., u n ) 
is invertible at 0. The inverse function theorem shows that U is a C 1 diffeo- 
morphism near 0. Moreover U is C 1+a and so is U _1 since 

DiU- 1 ) = (DU o t/" 1 )- 1 . (2.28) 

The fact that any C 1 diffeomorphism whose coordinate functions are A- 
harmonic is a C 1+a diffeomorphism follows from Proposition 12.61 □ 

Proof of Theorem \1.1[ As discussed above, in any local coordinate system 
where the metric has C a regularity, the p-harmonic equation is of the form 
div„4(x,V-u) = 0 where A is given by (12.21) and satisfies (12.31) (12.61) . The 
existence and regularity of p-harmonic coordinates is then a consequence of 
Theorem 12.11 for s < 1. The case s > 1 is contained in }LS14| . 

For s = 1, we have g £ W l ' q for all q < oo and thus by Proposition 12.71 
the p-harmonic coordinates are W 2,q regular. This regularity is enough for 
us to write the p-harmonic equation for each coordinate function u = u 3 in 
nondivergence form as 

a jk d jk u = f 

where the coefficients a? k and the function / are defined by 

a? k (x) = d^ k A j (x, Vu(i)), 
f(x) = -d Xj A J {x,Vu(x)). 

Now the linear equation a 3k dj k u = f is elliptic since Vtt is nonvanishing. 
Since u £ W 2 ' q for all q < oo, we have Vu £ C 7 for any 7 < 1 by the 
Sobolev embedding. Consequently a 3k £ C 7 for 7 < 1 and / £ C°. Thus 
u £ C *, see e.g. |Ta961 Theorem 14.4.2]. That the inverse of the p-harmonic 
coordinates are also of this regularity follows from the formula (|2.28l) . In de¬ 
riving this conclusion, note that Zygmund spaces are closed under pointwise 
multiplication |Tr83] . 

Suppose that the metric in original coordinates at xq is represented by the 
matrix Go, and let V be the corresponding A-harmonic coordinates. The 
metric in p-harmonic coordinates at xo is given by 

UP(0)" t G 0 T>P(0)” 1 . 

1 /2 

Choosing V so that DV( 0) is very close to Gq , we can arrange that the 
metric in p-harmonic coordinates is arbitrarily close to identity at xq. □ 

Corollary 2.8. Let ( M,g) be a Riemannian manifold with g £ W 1,q , q > 
n, in a local coordinate chart about some point xq £ M. There exists 
a p-harmonic coordinate system on a neighborhood of xo that is a local 
W 2 ' q diffeomorphism. Moreover, the Riemannian metric is W 1,q in these 
p-harmonic coordinates. 

Proof. Since W 1,q C C a for some a > 0, there exists a C 1+a regular p- 
harmonic coordinate system (v 1 ) on a neighborhood of xq by Theorem ll.il 
Let us denote V = (u 1 ,..., v n ). For g £ W 1,q , the function A defining 
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p-harmonic functions satisfies the conditions (12.201) and (12.211) . and thus by 
Proposition 12.71 we have that V £ W 2,q . By the formula 

DiV- 1 ) = (DV or 1 )' 1 

it is easy to see (by using integration by substitution and the standard 
cofactor formula for an inverse matrix) that the same holds for the inverse 
of V. Note that W 1,q is an algebra, see Section Q] 

That the Riemannian metric is W 1,q in p-harmonic coordinates follows 
from the coordinate transformation rule for the Riemannian metric 

~g = D(y- 1 ) T g\ v -.D(V- 1 ). 

Here g and g are the metrics in the original and in the p-harmonic coordi¬ 
nates respectively. □ 

3. Regularity of conformal and 1-quasiregular mappings 

The first application of n-harmonic coordinates is to the regularity of C 1 
conformal mappings, or W 1,n 1-quasiregular mappings, between Riemannian 
manifolds with C s , s > 0, s / 1, metric tensors. These results extend the 
corresponding results in [ Iw82. LSl4] that covered Holder exponents s > 1, 
and the earlier results [Le76[ IMa86] . 

The proofs are essentially the same as in [ LSI4 ] and are only sketched. We 
however cover a slight lapse in the proof of [ LS141 Theorem 3.1] concerning 
the integer values s = 2, 3,4,... of the regularity parameter s. We also 
comment on the case s = 1 after the theorem. 

Theorem 3.1. Let ( U,g ) and ( V,h ) be Riemannian manifolds such that 
g,h £ C s , s > 0, s / 1. If (f> : U —>■ V is a C 1 diffeomorphism that is 
conformal in the sense that <f*h = eg for some continuous positive function 
c, then 0 is a C s+l diffeomorphism from U onto V. 

Proof. Let first g,h £ C s with s > 0. Fix xq £ U, and let v be n-harmonic 
coordinates near 4>(xq). Since h £ C s , Theorem 11.11 ensures that v £ C* +1 . 
Let u = v o (f>. Since 4> is a C 1 conformal diffeomorphism, a straightforward 
calculation [LS14L Theorem 3.1] shows that u is an n-harmonic coordinate 
system near xo- Now g is C s , which implies that u is C'* +1 by Theorem ll.il 
Thus 0 = v~ 1 o u near xo and 0 -1 = it -1 o v near 4>{xq). 

Assume now that s > 0 is not an integer. Then u, v and their inverses 
are C s+1 , and thus f> and also f are C s+1 (the spaces C s+l = C* +1 are 
closed under composition for s > 0 noninteger). 

Assume next that s £ {2,3,4,...}. We first show that the conformal 
factor c in 

f>*h = eg 

is a C s function (this part of the proof is missing in ILS141 . Theorem 3.1]). 
The argument above yields that f £ C s+l ~ £ and c £ C s ~ e for any e > 0. 

Let us recall the transformation rules for the Ricci tensor and the scalar 
curvature under conformal scaling. If g = g, we have 

Rij(g) = Rij(g ) - (n - 2)(V i 3 j / - (■ dif^djf )) + (A/ - (n - 2)\df\l)g, 

R(g) = e~ 2f (R(g) + 2 (n - 1)A / - (n - 2)(n - l)\df\ 2 g ) (3.1) 
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if / and g are sufficiently smooth, see e.g. (Be87| . Since 4> is at least C 3 e 
and g , h are at least C 2 , we have 

g 2/ c — g~ 1 ( j)*h E C 2 ~ £ for any e > 0. 

This, together with the assumption that g, h E C 2 , yield that the identi¬ 
ties (13.11) hold in the distributional sense as described later in Section 0J In 
addition to the identities above, the tensoriality identities 

fRij(h) = Rijitfh), 

4>*R(h) = R(4>*h) 

also hold in the distributional sense. 

Let us now apply the scalar curvature to both sides of the equation 

4>*h = e 2 ^ g. 

The identities above yield 

4*R(h) = e~ 2 f(R(g ) + 2 (n - 1)A/ - (n - 2)(n - l)\df\ 2 g ). 

Since g,h E C s and 4> £ C s+l ~ £ , / E C s ~ £ , we may solve for A / which 
implies that A/ E C s ~ 2 . Now applying the Ricci tensor to (j)*h = e 2 ?g, 
using the identities above and using that A/ E C s ~ 2 , gives that 

dijf E C s ~ 2 for all i, j = 1,... n. (3.2) 

Thus /, and consequently c, is a C s function. 

The rest of the proof of the theorem now follows by expressing the second 
derivatives of 4> by the coordinate transformation formula 

didj<T = r *(c 5 )4<T - r]3(fc) 

of the Christoffel symbols (as in ('1170 . Ta()(i ). Since c E C s . the right hand 
side of this equation is in C s ~ l completing the proof. □ 

The reason we could not prove the theorem for s = 1 is that in this case 
the relation (G2D for the second derivatives of the conformal factor would 
involve terms containing second distributional derivatives of C 1 functions. 
We could not verify that the conformal factor is in C 1 in this case. 

The next result concerns the regularity of Riemannian 1-quasiregular 
mappings and the regularity of conformal bi-Lipschitz mappings. We re¬ 
fer to |DTT4l 111131 for details about Riemannian quasiregular mappings. 
The proof is the same as the proof of Theorem 4.4 in [ LSi4 ], except for the 
additional details for s = 2, 3,4, ... given above. We omit the proof. 

Theorem 3.2. Let ( M,g ) and ( N,h ) be Riemannian manifolds, n > 3, 
with g, h E C s , s > 0, s / 1. Let cf> : M —> N be a non-constant mapping. 
Then the following are equivalent: 

0 is a Riemannian 1-quasiregular mapping, 

4> is locally bi-Lipschitz and (f)*h = eg a.e., 

0 is a local C 1 diffeomorphism and (j)*h = eg, 

0 is a local C s+l diffeomorphism and (j)*h = eg. 
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Finally, we remark that the theorem holds also for n = 2 if we a priori 
assume to be a local homeomorphism. This gives an extension of I-ft 1 1 
Theorem 4.5] to Holder exponents s > 0, s / 1. 


4. WEYL TENSOR 


We now consider low regularity Riemannian metrics and establish an el¬ 
liptic regularity result for the Weyl tensor. We first assume that 

g jk € W 1 ’ 2 n L°° 

in some system of local coordinates near a point. (In this section all function 
spaces are assumed to be of the local variety near a point, that is, we write 
W 1 ' 2 instead of WK(U) etc.) This seems to be a minimal assumption for 
defining the Weyl tensor: the set W 1,2 fl L°° is an algebra under pointwise 
multiplication |KP 88 ], [ BBR12] , and therefore g jk , \g\ € W 1,2 n L°°. We see 
that 

r afc = \g cr (d a9 br + d b g ar - d r g ab ) € L 2 , 

Rabc = d a r d bc - d b r d ac + rixL - C c rL g w - 1 * 2 + l\ 

R bc = R abc a ew - 1 ’ 2 + L 1 . 

We write this schematically as 

Kb = Oig-'dg), 

R abc d = 0{d{g~ l dg ) + (< g~ l dgf ), 

R bc = 0{d{g~ l dg ) + (g^dg) 2 ) 

where 0(g~ 1 dg) denotes a finite linear combination of terms of the form 
g pq djg rs , etc. Now, since multiplication by an W 1,2 fl L°° function maps 
W 1,q to W 1,2 for any q > n, it also maps W~ 1,2 to W~ 1,q ', and we have 

Rabcd = O (gdig^dg) + g{g~ l dg) 2 ) 6 W~ 1,q ' + L 1 . 

Since 

Rab Tu 7T Rdab 

2{n — 1 ) 

where R = g rs R rs is the scalar curvature, we have 


Rab = 


n — 2 


abed — Rabcd T PacQbd Rhcfjad T Pbd 9 ac Podfjbc 

1 


(Rabcd T Rac 9 bd Rbcdad) T 
1 


n — 2 


1 ) (Racdbd RbcSad) 


+ 


' (Rbddac RadSbc ) 


R 


~(9ac 9 bd 9bc9ad) 


/ (n — i)(n — 2) 

= O (g{l + g~ l g) [d^dg) + (g~ l dg) 2 ]) € W~ 1,q ' + L 1 . (4.1) 


This shows that if gj k € I / F 1 ’ 2 nL°° in some system of local coordinates, then 
one can make sense of the components W abcc ] of the Weyl tensor as elements 
in W~ 1,q + L 1 in these coordinates for any q > n. It is easy to see that the 
identity W(cg) = cW(g) remains true if c, g ]k € W 1,2 fl L°°. 
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We wish to understand the equation W a b c d(g) = 0 as a quasilinear diver¬ 
gence form system for the matrix elements gjk- Below we write D = —iV. 

Lemma 4.1. Let gjk G W 1 , 2 nL 00 in some local coordinates. If W ab cd(g) = 0 
in these coordinates and if u G W 1 ’ 2 I ~1 L°° is the column vector consisting 
of the matrix elements gjk, then u is a distributional solution of 

Di(A lm (u)D m u) = B(u,Du) 

where B is given by 

B a (u, Du) = B l ™(u)Diu ■ D m u , 1 < a < M, 

and each A lm (resp. B l f n ) is an M x N (resp. N x N) matrix function whose 
entries are rational functions of the components of u. Moreover, each A lm 
and B 1 ™ is smooth in the set of vectors u corresponding to positive definite 
symmetric matrices. One can choose M = n 4 and N = n ( ra + 1 ) . 

Proof. Follows from (14.11) after using the Leibniz rule: 

W abcd = 0(d(g{ 1 + g~ 1 g)g~ 1 dg) - d(g( 1 + g~ l g))g~ x dg 
+ g(l+ g~ 1 g)(g~ 1 dg) 2 ). 

The fact that A lm and are rational functions follows by using Cramer’s 
rule for the inverse matrix (g^ k ). □ 

The main point, as observed in [LS15], is that the equation W a bcd{g ) = 0 
becomes elliptic if it is written in n-harmonic coordinates and the determi¬ 
nant is normalized to one. This is expressed by the next result. 

Lemma 4.2. Let gjk G W 1,g n C a in some local coordinates, where q > 2 
and 0 < a < 1. If g a b is the representation of g in any n-harmonic coordinate 
system, and if 

^lab = |S'| ^ gab, 

then g a b G W 1,q D C a . If u is the column vector consisting of the matrix 
elements g ab , and if A lm (x) = A i m 

(u(x)), then the linear operator 

V HA DiiA^DmV) 

is overdetermined elliptic (i.e. its principal symbol is injective). 

Proof. Since gjk G W l ' q P\C a , by Theorem 11.11 and Proposition 12.71 we know 
that any n-harmonic coordinate system is W 2 ' q D C 1+a . Now W 1,q D C a is 
an algebra under pointwise multiplication, which shows that g a b and g ab are 
in W 1,q C C a . 

To prove ellipticity, we need to show that the matrix valued symbol 
A lm £i£ m is injective for any £ G M ri \ {0}. In order to simplify notation 
we will omit all hats in the rest of this proof and write g a b instead of g ab etc. 
Assuming for the moment that g is C °°, we note the following consequences 
of Bianchi identities, 

V d M abc(i = 0, V a N abcd = 0 

where 

M a bcd — Babcd T Racgbd Phcfjad., Nabcd — Pafifjcd ^ 
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The left hand sides of these Bianchi type identities may be viewed as third 
order nonlinear operators acting on g. After linearization we obtain the 
following symbol versions of the previous identities, 

i d o<:(M abcd ) = 0, ti cl(T t( N abcd) = 0. (4.2) 

Here cr^(P) = a^(P)(g)h is the principal symbol of the linearization of an 
operator P , linearized at g = g a b and acting on some h = h ab . Now (E2D 
are algebraic equations for the principal symbols er ^{M abc d) and a^(N abc d) 
when g E C°°, but the since these principal symbols can also be computed 
for W l,q n C a metrics and their algebraic expressions remain identical, it 
follows that (14.211 is valid for g ab E W 1,q fl C a . 

By (14.111 we have 


^2 ^ 

^Pcibcd 44 abed V ^achd T 
n — 2 


n — 2 


(■Nu ac N a dbc ) 


n — 3 n — 3 

4 ~Pbc9ad 777 777 9ocf]hd T 


1 


n — 2 


2 (n — 2) 


2 (n — 2) 
R 


(n — l)(n — 2) 


{Pgbd9 ac RQadQbc) 

(dacQbd QbcQad )• 


Let us now assume that (J^{W a bcd)h = 0, where h = ( hij ). Contracting this 
equation by £ a £ d and using (|4.2|) yields 


Q = 2 (n = _ 2 |£|2 CT?(i4c)/i + 

O 77/ 77-1 

+ —■— 

n — 1 


The symbol on the right hand side is now (a scalar multiple of) the symbol 
of the Bach tensor, see [LS15, Lemma 3.2], The W l,q fl C a regularity is 
also enough for the n-harmonic coordinate gauge condition in }LS15| for 
the contracted Christoffel symbols (with g having determinant one) to be 
written as 

r k (9) = ^ g ka d a \ogg kk . 

The required ellipticity follows now from [LS151 Lemma 2.3]. □ 

Proof of Theorem \1.3[ The assumption that g E Ty 1,n nC' Q in some local co¬ 
ordinates together with Theorem 1 1.1 1 implies that there exists an n-harmonic 
coordinate system with C 1+a regularity. By Lemma 14.21 we know that in 
any n-harmonic coordinates, g ab and g ab = \g\~ 1 ^ n g a b are br W 1,n nC a . We 
will prove that 

9ab G C°°. (4.3) 

Since the Weyl tensor of g ab vanishes, we get that g ab = co5 a b where Co E C°°. 
Thus g ab = \g\ 1/n c 0 6 a b where c = \g\ l/n Co is a positive function in W 1 ’ n nC a 
as required. 

It remains to prove (14.311 . In the rest of the proof, we will omit hats and 
write g ab and u instead of g ab and u etc. Since the Weyl tensor of g vanishes, 
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Lemmas 14.1144.21 imply that u solves the overdetermined elliptic system 


D l (A lm (u)D m u) = B(u,Du) 


where u £ W 1,n n C a is the column vector consisting of the matrix elements 
of g a b. We claim that for any t > n, one has 


f D l (A lm (u)D m u) = B(u,Du) 

\ u£ w 1 ’* n c a 


U e w l ’ t+6 n C a 


(4.4) 


where (5 = If (14.41) is known, then starting from u £ W 1,n ft C a and 

iterating (14.41) gives u £ W 1,q for any q < oo. Then by Sobolev embedding 
u £ C 13 for any (3 < 1, and u solves 


D l {A lm {u)D m u) = B(u,Du) 


where A lm (u) £ C 13 and B(u,Du ) £ L q for all (3 < 1 and q < oo. The 
Schauder estimates in Proposition lA.II ci give that u £ C 1+/3 for all (3 < 1, 
and iterating the higher order Schauder estimates in Proposition lA.ll dl 
gives that u £ C°°. This shows (14.31) . 

To prove (14.41) . note that if u £ W 1,t nC a and t > n, then x i-a A lm (u(x)) 
and x i-A are in W l,t D C a (this set is an algebra under pointwise 

multiplication). Thus, if we write f(x) = B(u(x ), Du(x )), we see that u 
solves 

A(A ,m (u)D m u) = / 

with A lm (u ) £ W 1,t \~\C a and / £ L t ! 2 . The Calderon-Zygmund estimate in 
Proposition lA.ll bl implies that u £ IT 2,4 / 2 . But we also know that u £ C a , 
and C a C W £,q whenever e < a and q < oo. Under these conditions we 
have 

u £ w 2 ’ t/2 n w £ ’ q , 

and by complex interpolation u £ [VP 2 ’*/ 2 , W £,q ]g for 0 < 9 < 1. We choose 
6 so that 

(l-0)-2 + 0-e = l, 

which gives 0 = 2 ~ • Thus u £ W 1,p where 


1 

P 




i =n+,i. 

2 — e t q 


Since q < oo was arbitrary, we get u £ W 1,p whenever p < and e < a. 
In particular this gives (14.41) . □ 


Appendix A. Regularity for overdetermined elliptic systems 

Let C M n , n > 2, be a bounded open set. Let us consider an M x N 
linear system of PDEs P(x, D)u = f in R, where D = —iV and 

P{x,D)u = D l {A lm D m u) (A.l) 

and the coefficients satisfy 


A lm £ Ci oc (Q,R M xN ). 


(A.2) 
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Let p(x,£) = A lrn (x)^i^ m be the principal symbol of P. We assume that P 
is overdetermined elliptic in the sense that for any x G fl and £ G M n \ {0}, 
the matrix p(x,£) is injective: 

p(x,0 C = 0, (6R N =► C = 0. (A.3) 


If this holds, then necessarily M > N. 

We will prove the following C alder on- Zygmund and Schauder regularity 
results for overdetermined elliptic systems in divergence form. 


Proposition A.l. Let u G <f(a R' v ) solve P(x,D)u = / in the sense 
of distributions in f2, where P satisfies dSHD-dSUD and qo > 1. 

(a) If 1 < q < oo and 

A lm € c£ c , / G w£’ q , 

then u G W^ q - 

(b) If 1 < q < oo and 

A ,m € c» c n w&, f e L» oc , 

where 

J t = max{n, q}, q ^ n, 

\ t > n, q = n, 

then u G Ib^jL 

(c) If 0 < a < 1 and 

/^t—1+a 

^ u loc> / °Ioc > 

then u G C^ a . 

(d) If A: > 0 is an integer, 0 < a < 1, and 


a lm r- r-ik+cx. 
A e °loc > 


f€Ci 


k—1+a 
loc 


then u G Ct+ 1+a . 


Remark. As examples of functions in W^’ q and Cj o( | +Q , we mention that 
if 1 < q < oo, one has 



a < - JL - 

y ^ n—1 ’ 
^ > ri—T> 


a < 1, one has 


/ J e ifoc 




Remark. It is well known that for any q > 2 there is a uniformly elliptic 
scalar equation Dj(a qk D^u) = 0 where are L°°, such that there is a 
nontrivial VL 1,2 solution that is not in W l,q |Me63] . Thus (a) above does not 
hold if A lm G L°°, but (a) may remain true if one assumes A lm G L°°nVMO 
(see [DiF96] for the scalar case). 
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We will need some standard facts about function spaces, see [BL76j . [ Ta96l 
for more details. If s G R and 1 < p < oo, then W s,p = bF s,p (R n ) is the 
Banach space with norm 

ll/llir s .P(R") = \\{D) S f\\LP(R n )- 

Here, (£) = (1 + ICI 2 ) 1 ^ 2 and ( D) s is the Fourier multiplier corresponding to 
(£) s . In general, if m(£) is a C°° function which is polynomially bounded in 
R n together with its derivatives, the corresponding Fourier multiplier m{D ) 
is the map on 2/^(R n ) (the space of tempered distributions) defined by 

m(D)f = ^ -1 {m(0/(£)}- 

If k G {0,1, 2,...} and 0 < a < 1, then C k+a = C' fc+ “(R n ) is the Banach 
space with norm 

uni _ A"' ii 7-i7 #n i — D^f(y)\ 

\\f\\c k + a (R n ) — Y H-C )7 /IIl° 0 (R") + 2^ SU P-i—z—la-• 

| 7 |<fc Yi\=k x ^ y ' x y ' 

One has C k+a = C k+a = where Cl are the Zygmund spaces and 
the Besov spaces in R n . We write C~ l+a = C~ l+a for 0 < a < 1. The 
operator ( D )* is an isomorphism W S)P -A W s ~ t,p and Cl -A and Di 

is a bounded map W s,p —> W s ~ l ' p and Cl -A- C* -1 for any s,t G R and 
1 < p < oo [ Ta96l Section 13]. We will also use the Sobolev embeddings 

{ . L p *, 1 < p < n, 

OkocAoc’ P = ? b 

O 1 -^, n<p< oo, 

where p* = for 1 < p < n. 

The first step in the proof of Proposition IA.1I is a parametrix construction 
for solutions of P xo v = / in R n , where P xo the frozen coefficient operator 

P Xo = A lm (x 0 )Di m , xq G Q. 

We fix a cutoff function ^ G C£°(R n ) with tjj = 1 for |£| < 1/2 and -</ = 0 
for |£| > 1. 

Lemma A.2. Fix xo G fi. There is an operator E x0 = e Xo (D) such that for 
any v G S*"(M. n ,WL N ) one has 

v = E xo P xo v + ip(D)v (A.4) 

where e XQ : R n -A M. NxM is the bounded smooth matrix valued function 

e zo(0 = (! - V , (0)(K a; o,C)M a; o,0) _1 p( a; o,0*- 
The operator E Xo is a bounded map IF s,p (R n , R M ) -a- bF s+2,p (R n , l w ) and 
C*(R n ,R M ) —» C/ +2 (R n , R^) whenever s G R and 1 < p < oo. 

Proof. The operator P XQ has full symbol p(x o,£) = A lm (xo)£,i£ m . Then 

n 

p(x 0 ,0*p(x 0 ,0 = Y A:k (xo)*A lm (xo)^ k ^ m . 

j,k,l,m= 1 

The ellipticity condition (I A. 31) implies that 

det (p(x 0 , 0*P{ X o,0) 7^ °, £ G R n \ {0}. 
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Since p(xo; £)*p( x o> £) is homogeneous of degree 4 in £, the function e Xo 
defined in the statement of the lemma is smooth in £ and homogeneous of 
degree —2 for |£| large. 

If v € M. N ), the identity 

v(0 = (1 = e xo (£)p(x 0 , 0«(0 + 

implies v = E Xo P Xo v + ip(D)v by taking inverse Fourier transforms. Now 
e xo (£) is a matrix valued classical symbol in the class Sf 0 , and the mapping 
properties of the corresponding pseudodifferential operator [Ta96 l Section 
13] imply the mapping properties for E Xq . (Alternatively, one could apply 
the vector valued Mihlin multiplier theorem [BL761 Theorem 6.1.6] and its 
version on Besov spaces to rn XQ = (£) 2 e Xo .) □ 

Next we give a local representation formula for solutions of the equation 
P(x,D)u = f. To prove the Schauder estimates it is convenient to do a 
rescaling: if xq € 11 and r > 0 is small, we write 

A lm (y) = A lm (x 0 + ry), u(y) = u(x 0 + ry), f(y ) = f(x 0 + ry). 

We also fix a cutoff function x € C(? 0 (]R n ) with y = 1 for |x| <1/2 and 
supp(x) C £(0,3/4), and we let X 2 (y) = xiy/ 2). 

Lemma A.3. Let u € 1 / F 1 ) ) ' ( / 0 (1L M ,v ) solve P(x, D)u = / in Q, where go > 1- 
Fix a point xo € fl, and assume that £(xo,2r) C J). Then v = xu satisfies 
the equation 

v + E xo Q xo v = E X0 F + i/j(D)x 2 V in R n , 

where 

Q xo w = D l {xi{A lm - A lm (0))D m w) 

and 

F = r 2 X f + A lm D lX D m u + £,(A ,TO iZD m x). 

Proof. Note that u satisfies the equation 

D l {A lm D m u) = r 2 ~f in £(0,2). 

If v = x^i we compute 

DiiA^Dmv) = A(A* m £ m (x«)) = D^xDmU + A lm uD mX ) 

= X (r 2 f) + A lm D lX D m u + D l (A lm uD m x). 

Note that the left hand side operator frozen at 0 is just P XQ . Thus 

Di(A lm D m v) = P Xo v + Q Xo v 
where Q xo is as stated. It follows that 

P XQ v + Qx 0 v = F in M n . 


The result follows from Lemma I A. 21 


□ 
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Proof of Provosition \A. 1[ Let u G IF^’® (IL R iV ) solve P(x,D)u = f in 0, 
where qo > 1. Fix xo G 12, and assume that B(x q, 2 r) C 12. By Lemma fA.31 
the function v = satisfies 

T(v) = G in R n , 

where 

T(w) = w + E Xo Q Xo w, 

G = E X0 F + ij}{D)x 2 V. 

Assume now the conditions in (a). We may assume qo < n, so IF 1,90 C L q o 
by Sobolev embedding. We may also assume qo < Q (otherwise there is 
nothing to prove). Fix po > 2 so that qo,q G \p' 0 ,po]- We claim that if r is 
chosen small enough, then 

j T is invertible on IF 1,p (R n , M N ) for p G [p' 0 ,po\, and (A.5) 

\ G G W 1,p for p G [q 0 , minj^Q, g}]. (A.6) 

Assuming these two claims, it follows that v = T~ 1 (G) G W 1,rain ^ q 0 ’ 9 i. If 
g < gj we have v G W 1,q and we are done. Otherwise we get v G IF 1,9 o, 
and repeating the argument with qo replaced by q$ we see that G and thus 
also v are in WA> min {?o ,9 i. Iterating this procedure finitely many times gives 
v G IF 1,9 , which proves u G W^ q upon varying xq. To see (1A.5D . note that 

\\E xo w\\ w i,p < C\\w\\ w -i, P , w G IF _1 ’ P , p G \p' 0 ,Po\- 

This follows by complex interpolation from the cases p = po and p = p' 0 
given by Lemma I A. 21 Consequently 

\\E Xo Q Xo w\\ wl , P < C\\Di( X 2 (A lm - A lm (0))D m w)\\ w -r P 

1,771 

Since A lm are continuous, we can choose r small enough so that 

\\E X0 Qx 0 w\\wrp < \\\w\\w^i P € \p'o,Po\- 

This implies (|A.5D by Neumann series. To see (1A.6D . note that our assump¬ 
tions and Sobolev embedding imply 

F G IF -1 ’ 9 + L qo + IF _1 ’ 9 °*. 

Since all the terms are compactly supported, we get 

F G W~ 1,p , p G (1, min{gg, g}]. 

The mapping properties of E Xq in Lemma lA.2l and the fact that ip(D)x 2 ^ is 
Schwartz imply (IA.6I1 . 

Assume next the conditions in (b). We claim that for r > 0 small, 
j T is invertible on IF 2 ’ 9 (R n , R' v ), and 
1 G G IF 2 ’ 9 . 


(A.7) 

(A.8) 
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Since also the conditions in (a) hold, we have v = T~ l (G) and thus (b) 
follows immediately from these claims upon varying x$. To prove (IA.7p . 
note that Lemma lA.21 gives 

\\ E xo w \\w 2 ’i < C\\w\\ Lq , w £ L q . 

Combining this with the pointwise multiplier property 

||au||vyi, 9 < ||art||i,9 + || V(<m)||z,? < C(||a||Loo + ||Va|| L t)||rt|| w i, g 

by using the assumption on t (if p = n this holds for a supported in H(0, 2)), 
we get 

n 

\\E X0 Q X0 w\\ W 2 , q < C]T||x 2 (i Zm - A lm (0))D m w\\ w i, q 
1=1 

< Co r ^Q(l)\\w\\ W 2, q 

using the fact that A lm € C® oc D This shows (1A.7I) by Neumann series. 
To prove (IA.8I) . consider first the case 1 < q < n. Then L q C W~ 1,q , 
so u € Wfof! by (a). Since L n ,L q * C L q , we get F € L q and therefore 
G € W 2,q by the mapping properties of E Xq . In the case n < q < oo we 
get L q oc C Wj~p’ p for all p < oo, thus u £ W\ 0 <‘ for all p < oo by (a), and 
u £ L°° by Sobolev embedding. The assumption A lm £ C® oc D gives 
F £ L q and G £ W 2,q . This proves (|A.8D and also (b). 

Assume next the conditions in (c). It is enough to prove that for r > 0 
small, 

(T is invertible on C 1+Q (M n ,R^), and (A.9) 

1 G £ C 1+a . (A.10) 

The main point is that \A lm (y) — A lm (z)\ < M\y — z\ a , so we have 

\\ X2 (A lm -A lm m\c‘><Cr a 


and consequently 


\\E Xo Q Xo w\\c^ — C 

i=i 


< C'r a ||u!|| (: 7i+ a 


A lm (0))D m w\\c ■ 


This proves (IA.9D by Neumann series if r > 0 is small enough. To prove 
(IA.10D one first notes that / £ ,q for all q, hence u £ W^ Q q for all q by 
(a) and u £ C^ oc for 0 < /? < 1 by Sobolev embedding. Thus F £ C~ l+a , so 
that G £ C 1+a by the mapping properties of E xo . 

Finally, assume the conditions in (d). We do the proof by induction on 
k. If k = 0 the result follows from (c). Assume the result holds smoothness 
indices < k — 1 where k > 1, and assume that 


Aim /~ik+a 

A ^Toc ’ 


/eq 


k— 1+a 
loc 


By induction hypothesis, u £ for f3 < 1. Also, since k > 1, part (b) 

implies that u £ W^ q for any q < oo. Hence it is possible to differentiate 
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the equation once to get 

D^DmiDju)) = Djf - DtdDjA^Dmu). 

The right hand side is in C^~ 2+a , hence the induction hypothesis gives 
DjU £ Cf 0 + Q for all j. This proves that u £ C^ 1+a . □ 
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